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Abstract. We consider diffusion processes xt on the unit interval. Doob- 

transformation techniques consist of a selection of xj— paths procedure. The 

law of the transformed process is the one of a branching diffusion system 

of particles, each diffusing like a new process xt, superposing an additional 

drift to the one of xt- Killing and/or branching of x^— particles occur at 

^ y ' some space-dependent rate A. For this transformed process, so in the class of 

• ^ branching diffusions, the question arises as to whether the particle system is 

O 1 sub-critical, critical or super-critical. In the first two cases, extinction occurs 

with probability one. 

We apply this circle of ideas to diffusion processes arising in population 
^^^ genetics. In this setup, the process xt is a Wright-Fisher (WF) diffusion, 

either neutral or with mutations. 

We study a particular Doob transform which is based on the exponential 
function in the usual fitness parameter a. We have in mind that this is an 
1^ alternative way to introduce selection or fitness in both WF-like diffusions, 

N] , leading to branching diffusion models ideas. For this Doob-transform model 

>«.i^ ' of fitness, the usual selection drift ax (1 — x) should be superposed to the one 

VJ ' of xt to form xt which is the process that can branch, binarily. 

CN In the first neutral case, there is a trade-off between branching events giving 

^^^ , birth to new particles and absorption at the boundaries, killing the particles. 

^"^ ' Under our assumptions, the branching diffusion process gets eventually glob- 

ally extinct in finite time with exponential tails. 

In the second case with mutations, there is a trade-off between killing events 
removing some particles from the system and reflection at the boundaries 
where the particles survive. This branching diffusion process also gets eventu- 

1 ^ ■ ally globally extinct but in very long finite time with power-law tails. 
K> ' Our approach relies on the spectral expansion of the transition probability 

3 , kernels of both xt and xt- 



r^ 
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1. Introduction 

We consider diffusion processes on the unit interval with in mind a series of elemen- 
tary stochastic models arising chiefly in population dynamics. Special emphasis is 
put on Doob-transformation techniques of the diffusion processes under concern. 
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2 THIERRY E. HUILLET 

Most of the manuscript's content focuses on the specific Wright-Fisher (WF) dif- 
fusion model and some of its variations, describing the evolution of one two-locus 
colony undergoing random mating, possibly under the additional actions of muta- 
tion and selection. These models found their way over the last sixty years, chiefly 
in mathematical population genetics. We refer to the general monographs [3], |21) . 

[g, [g and [g. 

We now describe the content of this work in some more details. 

Section 2 is devoted to generalities on one-dimensional diffusions on the unit interval 
[0,1], say {xt'it > 0). Special emphasis is put on the Kolmogorov backward and 
forward equations, while stressing the crucial role played by the boundaries in such 
one-dimensional diffusion problems. Some questions such as the meaning of speed 
and scale functions, existence of an invariant measure, random time change... are 
addressed in the light of the Feller classification of boundaries. When the boundaries 
are absorbing, the important problem of evaluating additive functional a along 
sample paths is then briefiy discussed, emphasizing the prominent role played by 
the Green function of the model. 

So far, we have dealt with a given diffusion process Xt, and recalled the various 
ingredients for computing the expectations of various quantities of interest, sum- 
ming up over the history of its paths. In this setup, there is no distinction among 
paths with different destinations, nor did we allow for annihilation or creation of 
paths inside the domain. The Doob transform of paths is an invitation to do so. 
This important class of transformations is a particular instance of a more general 
construction based on multiplicative functionals. We fix the background. 

Roughly speaking, the Doob transformation of paths procedure allows to select 
sample paths x ^ y within any laps of time i, favoring large values of the ratio 
a (y) /a {x) for some specific functional a > that fixes the selection problem 
under study. The process solving this selection of paths procedure belongs to a 
class of branching diffusion processes, where independent particles diffusing like a 
new process Xt inside the interval are allowed to duplicate would the visited region 
of the state-space fulfill the selection of paths criterion or to die, if not. In the 
process, advantageous regions of the state-space are reinforced while unfavorable 
ones are left unexplored which is a reasonable physical way to look at selection of 
paths. The new process Xt alluded to is obtained from Xt just after superposing an 
additional suitable drift to the latter process. An important parameter is the state- 
dependent rate A at which killing and/or branching occur. Depending on A and on 
the type of boundaries which {0, 1} are to a;*, the full transformed process can have 
two stopping times: the time to absorption at the boundaries and the killing time 
inside the domain. Besides, there is or not an opportunity that the Xt— particles 
duplicate, leading or not to a daughters particle system evolving independently 
starting from where their mother particle died. The killing/branching issues depend 
on the sign of A. 

It turns out that the same diffusion methods used in the previous discussion on 
simple diffusions apply to the transformed processes obtained after the induced 
change of measure. We develop this circle of ideas. 
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We next apply these general ideas to diffusion processes arising in population ge- 
netics. 

In Section 3 we start recalling that WF diffusion models with various drifts are 
continuous space-time models which can be obtained as scaling limits of a biased 
discrete Galton- Watson model with a conservative number of offsprings over the 
generations. Sections 4 and 5 are devoted to a detailed study of both the neutral 
WF diffusion process (WFN), the WF diffusion with selection (WFS), the WF dif- 
fusion with mutations (WFM) and the WF diffusion with mutations and selection 
(WFMS) respectively. 

In this context, our suggestion is the following one: we can view the action of 
selection (or fitness) on the evolution of the allele frequency distribution, either 
neutral or with mutations, as a functional deformation of the sample paths of the 
original process, say Xt-, favoring initial values xq = x with small a (x) and terminal 
values Xt — y with large a (y) , for each t. In our construction, a (x) — e"'^ , cr > 0, 
is the chosen exponential fitness functional. Stated differently and more precisely, 
ii p{x;t,y) is the transition probability density of xt (either WFN or WFM), our 
model of the action of fitness is 

fitness 6°'*' 

p{x;t,y) -^'"'p{x;t,y)^-—p{x;t,y). 

With this choice of a, the modification consists of selecting those paths x ^ y oi 
Xt for which e'"^^"^' is large. As a result of this transformation of paths, the usual 
positive selection drift ax (1 — x) has to be superposed to the one of xt to form the 
new process Xt , but our functional definition of fitness also generates an additional 
branching multiplicative term, translating that a particle system pops in: The 
resulting transformed process is a (binary) branching diffusion of WF diffusions Xt ■ 
We may call the obtained processes the branching neutral Wright-Fisher process 
and the branching Wright-Fisher process with mutations. This point of view seems 
to be new, to the best of the author's knowledge. 

Because the spectral representation of both transition probability densities of WFN 
or WFM are known explicitly from the works of Crow and Kimura (see [U] , [3] and 
[3]), some easy consequences on the spectral structures of the branching transformed 
processes are available. For instance, it is possible to decide whether the BD process 
is sub-critical, critical or super-critical in the sense of ([1] and [2]). 

In Section 6, we therefore give a detailed study of the binary branching diffusion 
process obtained while using the Doob transform a (x) = e'^^ when the starting 
point process is a WFN diffusion process. We end up with a branching particle sys- 
tem, each diffusing according to the WF model with a selection drift, but branching 
at a bounded rate 6 > 0. In this setup, the particles cannot get killed, rather they 
are allowed either to survive or to split: the transformed process is a pure binary 
branching diffusion. For this super-critical binary branching diffusion process, there 
is a trade-off between branching events giving birth to new particles and absorp- 
tion at the boundaries, killing the particles. Thanks to the spectral representation 
of the WFN process, this problem is amenable to the results obtained in ([T] and 
[2]). Under our assumptions, the branching diffusion process turns out to be glob- 
ally sub-critical: the branching diffusion process gets eventually globally extinct in 
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finite exponential time. This requires the computation of the ground states associ- 
ated with the smallest nonnegative eigenvalue of the infinitesimal generator of the 
transformed process which are here shown to be explicit. In particular, the expres- 
sion of the quasi-stationary distribution of the particle system can be obtained in 
closed-form. 

In Section 7, we study the binary branching diffusion process obtained while using 
the same Doob transform, when the starting point process is now a WF diffusion 
process with mutations, assuming reflecting boundaries. We end up in a branching 
particle system, each diffusing according to the WF model with a mutation and 
selection drift, but branching at quadratic rate A, which is bounded from below 
and above. Although the particles are still allowed to split, they can now also get 
killed at the branching times: the transformed process is again a binary branching 
diffusion but with killing now allowed. In this setup, there is a competition between 
branching/ killing events and reflection at the boundaries where the particles sur- 
vive. This problem is also amenable to the results obtained in ([T] and [5]) and we 
end up now in a globally critical branching particle system, each diffusing accord- 
ing to the WF model with a mutation and selection drift. This branching diffusion 
process turns out to be globally critical: it also gets eventually globally extinct but 
now in long finite time, with power-law tails. 

2. Diffusion processes on the unit interval and Doob transforms 

We start with generalities on one-dimensional diffusions with the WF model and 
its relatives in mind. For more technical details, we refer to [5], [7], [13] and |19) . 
We also introduce Doob transforms as particular instances of the modification of 
the original diffusion process through a multiplicative functional. 

2.1. One-dimensional diffusions on the interval [0,1]. Let {wt;t>0) be a 
standard one-dimensional Brownian (Wiener) motion. We consider a 1— dimensional 
Ito diffusion driven by {wt;t > 0) on the interval say [0, 1], see [11]. We assume it 
has locally Lipschitz continuous drift / (x) and local standard deviation (volatility) 
g (x), namely we consider the stochastic differential equation (SDE): 

(1) dxt = / (xt) dt + g (xt) dwt, xq ^ x e I :^ (0, 1) . 

The condition on / (x) and g (x) guarantees in particular that there is no point x* 
in / for which \f {x)\ or \g {x)\ would blow up and diverge as [a; — a;, | — > 0. 

The Kolmogorov backward infinitesimal generator of ([T]) is G = / {x) dx + ^g"^ {x) d^. 
As a result, for all suitable ip in the domain of the operator St := e**^, u := u {x, t) = 
Eip {xtf\T^ ) satisfies the Kolmogorov backward equation (KBE) 

dtu = G {u) ; u{x,0) — ip (x) ■ 

In the definition of the mathematical expectation m, we have t At^ := inf (t, r^;) 
where Tx indicates a random time at which the process should possibly be stopped 
(absorbed), given the process was started in x. The description of this (adapted) 
absorption time is governed by the type of boundaries which dl :— {0, 1} are to 
{xt;t>0). A classification of the boundaries exists, due to Feller (see [13] pp. 
226): they can be either accessible (namely exit or regular), or inaccessible (namely 
entrance or natural). 
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2.2. Natural coordinate, scale, speed measure, time change. For such Mar- 
kovian diffusions, it is interesting to consider the G— harmonic coordinate (p G C^ 
belonging to the kernel of G, i.e. satisfying G {(p) — 0. For tp and its derivative 
(fi' := dip/dy, with (xq, j/o) G (0, 1), one finds 



f'iy) = ¥''(j/o)e ■'yo-^^M 

fix) = LP (xn) + ip' (yo) e ''■'^^^(^"^(iy. 



-2 P -^d^ 

J Xo 

One should choose a version of (/3 satisfying ip' (y) > 0, y E I. The function p 
kills the drift / of {xt;t > 0) in the sense that, considering the change of variable 
yt = p{xt), 

dyt = (p'g) {p^^ (yt)) dwt, yo^/p (x) . 
The drift-less diffusion (y*; i > 0) is often termed the diffusion in natural coordinates 
with state-space [ip> (0) , p (1)] =: p (/). Its volatility is g {y) := {(p'o) {f'^ iv)) ■ The 
function p is often called the scale function. 

Whenever p (0) > — cxd and (p (1) < +oo, one can choose the integration constants 
defining ip (x) so that 

rx -2fy4^dz, 

y^(^)= ,1 ^2fy4^dz ' 

with ((5(0) = and (^3(1) = 1. In this case, the state-space of {ytlt > 0) is again 
[0, 1] , the same as for {xt'.t > 0) . 

Finally, considering the random time change t -^ Ot with inverse: 9 -^ tg defined 
by Oti, — and 

e= f'<f{ys)ds, 
Jo 
the novel diffusion (wg := ytg',0 > 0) is easily checked to be identical in law to a 
standard Brownian motion on p) (/) . The random time tg can be expressed as 

tg = m (p^'^ (wr)) {(p^^y {wr) dr 
Jo 

where m{x) := 1/ yg^ip'j {x) is the (positive) speed density at a; = p^^ (y). Both 

the scale function p and the speed measure dfx = m (x) ■ dx are therefore essential 

ingredients to reduce the original stochastic process {xt;t>Q) to the standard 

Brownian motion {wg ;9 > 0) . The Kolmogorov backward infinitesimal generator 

G may then be written in Feller form 



2 dfj, \dip 



Examples (from population genetics): 

• Assume f (x) = and g'^ (x) = x {1 — x). This is the neutral WF model discussed 
at length later. This diffusion is already in natural scale and p (x) — x, m (x) = 
[x (1 — x)]~ . The speed measure is not integrable. 

• With TTi, 7r2 > 0, assume / (x) — tti — (tti -1- 7^2) x and g^ (x) — x(l — x). This 
is the WF model with mutation. The parameters 7ri,7r2 can be interpreted as 
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mutation rates. The drift vanishes when x — tti/tt (where vr := tti + 7r2 is the total 
nnitation pressure) which is an attracting point for the dynamics. Here: 

^' (y) - v' (2/0) 2/-^-^ (1 - y)~'"^ .v{x)^v M + v' {yo) ^ v'^^' (1 - y)~'^' dy, 

with (fi (0) = — cx) and ip (1) = +00 if 7ri,7r2 > 1/2. The speed measure density is 
m (x) oc x^'^i"^ (1 — x) ^'^~ and so is always integrable. After normalization to 1, 
m{x) is the beta(27ri, 2712) density. 

• With (7 G R, assume a model with quadratic logistic drift / (a;) = ax{l — x) and 
local variance g^ {x) = x{l — x). This is the WF model with selection or selection. 

For this diffusion (see [l^), if [x) = \n'e-^'^ ^'^'^ ™ (^) °*- [^ (1 ~ ^)]~ ^^'^'^ i^ ^^^ 
integrable. Here, u is a selection or fitness parameter. 

Time change and subordination. We start from the diffusion (1) with infini- 
tesimal generator G = fd^ + \g^&l and consider the time change problem without 
passing first in natural coordinate. Let the random time change 



t-^9t^ / g^{xs)ds 
Jo 





to 



Its inverse: 8 -^ tg defined by 9tg = 6* is given by 6* = L" g^ (x^) ds. 

In this new stochastic time clock, the subordinated diffusion {yg := xt^ ; > 0) obeys 

FF)' 



the Langevin SDK with potential U (y) := -2 J^ ■^M)'^^ 



f 
dye = ^ {yg)d9 + dwe, 
9 

with backward infinitesimal generator G — g^^G — -i^dx + ^9| (See 
164-169). 



pp. 



We have 9t — g^ (xt) meaning that at each point xt of the former motion, the 
motion of the path is accelerated or decelerated, depending on the rate g"^ {xt) ^ 1. 

Note that conversely tg = 1/g^ {yg) . Under the time substitutions, the road maps 
of the paths of both {xt]t > 0) and {yg;9 > 0) remain exactly the same. If a path of 
the former process is accelerated or decelerated by its squared volatility g^ (its local 
variance) at each locality, then this process boils down to the latter one. Stated 
differently, if we measure time by the amount of squared volatility accumulated 
within each of its path, the process {xt;t > 0) becomes {yg;9 > 0), both with state- 
space /. 

2.3. The transition probability density. Assume that / (x) and g (x) are now 

differentiable in /. Let then p{x;t,y) stand for the transition probability density 
function of Xt at y given xq = x. Then p := p {x; t, y) is the smallest solution to the 
Kolmogorov forward (Fokker-Planck) equation (KFE) : 

(2) dtp^G*{p),p{x-Q,y)^5y{x) 

where G* (•) ^ -dy (/ (y) •) + ^d^ {g^ (y) •) is the adjoint of G (G* acts on the 
terminal value y whereas G acts on the initial value x) . The way one can view this 
partial differential equation (PDE) depends on the type of boundaries that {0, 1} 
are. 
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Suppose for example that the boundaries o := or 1 are both exit (or absorbing) 
boundaries. From the Feller classification of boundaries, this will be the case if 
Vyoe(0,l): 

(3) (i) m{y) ^ Li{yo,o) and (m) ip' (y) / m{z) dz e Li{yo,o) , 

Jya 

where a function / (y) G ii (j/o, °) if /J|, 1/ {y)\ dy < +oo. 

In this case, a sample path of {xt ]t > 0) can reach o from the inside of / in finite time 
but cannot reenter. The sample paths are absorbed at o. There is an absorption 
at o at time t^^o = inf {t > : Xt = o \ xq = x) and P {tx^o < oo) = 1. Whenever 
both boundaries {0, 1} are absorbing, the diffusion xt should be stopped at t^ '■= 
Tx,Q Ar^j^i. Would none of the boundaries {0, 1} be absorbing, then t^ = +oo. This 
occurs when the boundaries are inaccessible. 

Examples of diffusion with exit boundaries are the WF model and the WF model 
with selection. In the WF model including mutations, the boundaries are entrance 
boundaries and so are not absorbing. 

When the boundaries are absorbing, then p (a;; i, y) is a sub-probability. Letting 
Pf{x) :— Jq p{x;t,y) dy, we clearly have Pt{x) — P {t^ >t). Such models are 
non-conservative. 

For one-dimensional diffusions, the transition density p (x; i, y) is reversible with 
respect to the speed density (|13j. Chapter 15, Section 13) and so detailed balance 
holds: 

(4) m (x) p {x; t,y) ^m (y) p {y;t,x) , < x,y < 1. 

The speed density m (y) satisfies G* (m) = 0. It may be written as a Gibbs measure 
with density: m{y) oc -j] ^ e~^^^^ where the potential function U (y) reads: 



(5) 



U{y) = -2 r IS^dz,0<y<l 
Jo 9 [z) 



and with the measure -^tt-\ standing for the reference measure. 

Furthermore, if p (s, x; t, y) is the transition probability density from (s, x) to (i, y), 
s < t, then —dsP — G{p), with terminal condition p(t,x;t,y) = Sy{x) and so 
p (s, x; t, y) also satisfies the KBE when looking at it backward in time. The Feller 
evolution semigroup being time-homogeneous, one may as well observe that with 
p := p (x; t, y), operating the time substitution t — s -^ t, p itself solves the KBE 

(6) dtp = G{p), p{x;0,y) = Sy{x). 

In particular, integrating over y, dtPt {x) — G (pj (x)), with p^ (x) = 1 (x G (0, 1)). 

p (x; i, y) being a sub-probability, we may define the normalized conditional proba- 
bility density q (x; i, y) :— p (x; t, y) / p^ (x), now with total mass 1. We get 

dtq = -dtPt {x) /pt {x)-q + G* (q) , q (x; 0, y) = 6y (x) . 

The term bt (x) := —dtPt (x) j p^ (x) > is the time-dependent birth rate at which 
mass should be created to compensate the loss of mass of the original process due 
to absorption of {xt\t > 0) at the boundaries. In this creation of mass process, a 
diffusing particle started in x dies at rate bt (x) at point (t, y) where it is duplicated 
in two new independent particles both started at y (resulting in a global birth) 
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evolving in the same diffusive way u. The birth rate function bt {x) depends here 
on X and t, not on y. 

When the boundaries of Xt are absorbing, the spectra of both —G and —G* are dis- 
crete (see [13] pp. 330): There exist positive eigenvalues {Xk)k>i ordered in ascend- 
ing sizes and eigenvectors (wfe,Mfe)fc>i of both —G* and — G satisfying —G* (vk) = 
XkVk and ~G {yk) = \kUk such that, with {uk,Vk) := ^^ Uk{x)vk{x)dx and 
hk '.= {uk,Vk) , the spectral representation 

(7) p{x;t,y) = Y,bke-^''*Ukix)vk{y) 

k>l 

holds. 

Let Ai > Ao = be the smallest non-null eigenvalue of the infinitesimal generator 
— G* (and of — G). Clearly, — jlogpj (x) — > Ai and by L' Hospital rule therefore 

bt (x) — > Ai. Putting dtq = in the latter evolution equation, independently of 

the initial condition x 

(8) q{x;t,y) ~> q^ (y) ^ vi (y) , 

where vi is the eigenvector of — G* associated to Ai, satisfying — G*wi — AiWi. 
The limiting probability wi/norm (after a proper normalization) is called the quasi- 
stationary Yaglom limit law of {xt ]t > 0) conditioned on being currently alive at 
all time t (see [23]). 

2.4. Additive functionals along sample paths. Let {xt;t > 0) be the diffusion 
model defined by ([T]) on the interval / where both endpoints are assumed absorbing 
(exit). This process is thus transient and non-conservative. We wish to evaluate 
the nonnegative additive quantities 



a (x) = E I / c{xs) ds + d{xr-,) 



where the functions c and d are both assumed nonnegative on / and dl — {0, 1}. 
The functional a (x) > solves the Dirichlet problem: 

-G(a) = c if a; £ / 
a ~ d if X G 9/, 

and a is a super-harmonic function for G, satisfying — G {a) > 0. 

Some examples: 

1. Assume c = 1 and d = : here, a — E(tj^) is the mean time of absorption 
(average time spent in (0, 1) before absorption). 



Consider a diffusion process with forward infinitesimal generator G* governing the evolution 
of p{x; t, y) . Suppose that a sample path of this process has some probability that it will be killed 
or create a new copy of itself, and that the killing and birth rates d and b depend on the current 
location y of the path. Then the process with the birth and death opportunities of a path has the 
infinitesimal generator A (j/) • +G* (■) , where A (y) = b (y) — d (y). The rate can also depend on t 
and X. 
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2. Whenever both {0, 1} arc exit boundaries, it is of interest to evahiate the 
probabihty that xt first hits [0, 1] (say) at 1, given xq — x. This can be obtained 
by choosing c — Q and d (o) = 1 (o = 1) . 

3. Let y & I and put c = ^1 (x e (y — e, j/ + e)) and d = 0. As e — > 0, c con- 
verges weakly to Sy{x) and, a =: g (x, y) — E (hm ^ /^ ^ l(j,_g j^^^) (x^) ds) = 
L p (x; s, y) ds is the Green function, solution to: 

-G (o) = Sy (x) ii X e I 
= Oiixedl. 

Q is therefore the mathematical expectation of the local time at y, starting from x 
(the sojourn time density at y). The solution is known to be (see |13j . pp. 198 or 
0, pp. 280) 

Q{x,y) = 2m{y) (^(1) -<^(y)) '^ ^.^j ' ^ {"j if x < y 

(y9(l)~(y9(0) 

(9) (x, y) = 2m (y) (^ (y) - y, (0)) '^^,|^~^^,g if x > y. 

The Green function is of particular interest to solve the general problem of evalu- 
ating additive functionals a (x) . Indeed, as is well-known, see [13] for example, the 
integral operator with respect to the Green kernel inverts the second order operator 
—G leading to 



a(x) = g{x,y)c{y)dy iix el 

a — d if X e dl. 

Under this form, a (x) appears as a potential function and any potential function is 
super-harmonic. Note that for all harmonic function h > satisfying —G (h) — 0, 



ah (x) :== / (x, y) c (y) dy + h (x) 
is again super-harmonic because — G (ah) = c > 0. 

2.5. Transformation of sample paths (Doob-transform) producing killing 
and/or branching. In the preceding Subsections, we have dealt with a given pro- 
cess and recalled the various ingredients for the expectations of various quantities 
of interest, summing over the history of paths. In this setup, there is no distinction 
among paths with different destinations nor did we allow for annihilation or cre- 
ation of paths inside the domain before the process reached one of the boundaries. 
The Doob transform of paths allows to do so. 

Consider a one-dimensional diffusion (xt; i > 0) as in ([T]). Let p (x; i, y) be its tran- 
sition probability. Let a (x) > as x G [0, 1] . 

Define a new transformed stochastic process (xf ]t > 0) hy its transition probability 

(10) p{x;t,y) ^ —--p{x;t,y). 

a (x) 
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In this construction of {Tt;t > 0) through a change of measure, sample paths 
X -^ y oi [xt] t > 0) with a large value of the ratio a (y) /a (x) are favored. This is 
a selection of paths procedure due to Doob (see [B]). 

The KFE for p clearly is dtp ~ G (p), with p{x;0,y) = 5y[x) and G (p) = 
a (y) G* {p / a (y)) . The adjoint Kolmogorov backward operator of the transformed 
process is therefore by duality 

(11) G{-) = ^Giaix)-). 

a (X) 

Developing, with a' (x) :— da (x) /dx and G (•) := ^g^dx (•) + G (•), we get 

(12) G(-) = -G{a) ■ +G{-) =: A(x) • +G {■) 

a 

and the new KB operator can be obtained from the latter by adding a drift term 
—g^dx to the one in G of the original process to form a new process {xt\t > 0) 
with the KB operator G and by killing or branching its sample paths at rate 
A (x) := G (a) /a. In others words, with / (x) := f (x) + —g^ (x) , the novel time- 
homogeneous SDE to consider is 

(13) dxt = / (xt) dt + g (xt) dwt, xq ^^^ x e (0, 1) , 

possibly killed or branching at rate A (x) as soon as A 7^ 0. Whenever {xt;t > 0) 
is killed, it enters conventionally into some coffin state {d} added to the state-space. 

Let us look at special cases: 

(i) Suppose a > is such that —G (a) > (By a > 0, we mean a > in /, possibly 
with a (0) or a(l) equal 0). Then a is called a super-harmonic (or excessive) 
function for the process with infinitesimal generator G. 

In this case, the rate X{x) —: —d{x) satisfies \{x) < and only killing occurs 
at rate d{x). Let Tx be the new absorption time at the boundaries of {xt\t > 0) 
started at x (with Tx — 00 would the boundaries be inaccessible to the new process 
Xt). Let Tx,d be the killing time of {xt]t > 0) started at x (the hitting time of 
d), with Tx,d = 00 if G (a) = 0. Then Tx ■= Tx A Tx.d is the novel stopping time 
for {xt;t > 0) . The SDE for {xt]t > 0), together with its global stopping time Tx 
characterize the new process {xt;t > 0) with full generator G to consider. 

(a) Suppose a > is such that — G (a) < 0. Then a is called a sub-harmonic 
function for the process with generator G. 

In this case, the rate A (x) =: b (x) satisfies A (x) > and only branching occurs 
at rate b (x) . The transformed process (with infinitesimal backward generator G) 
accounts for a branching diffusion where a diffusing mother particle (with generator 
G and started at x) lives a random exponential time with constant rate 1. When 
the mother particle dies, it gives birth to a spatially dependent random number 
M (x) of particles, with mean fJ,{x) = 1 + A (x) (where M (x) = 1 + A(A(x)) 
and A(A(a;)) is a geometrically distributed random variable on {0,1,2,...} with 
mean A (x)). Then M (x) independent daughter particles are started afresh where 
their mother particle died, with the event M (x) = impossible; they move along a 
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diffusion governed by G and reproduce, independently and so on for the subsequent 
particles. 

[iii) If A (x) =: b {x) is bounded above, A may be put under the alternative form 

A(x) = A* {p(x) - 1) , 

where A* = sup^grp ij A (cc) and 1 < /i (a;) < 2. In this case, we can assume that 
M (x) can only take the values 1 or 2 with probability pi (x) and p2 {x) respectively, 
with pi (x) +p2{x) ^l. Then, ^ (x) == E (M (x)) = pi (x) + 2p2 (x) = 1 + p2 (x) 
and 

A (x) = A*p2 {x) . 
Note that /i (x) — 1 — P2 {x) > Pi (x) = 2 — /i (x) . We get a binary branching process 
at rate A» with the event to produce two particles being more likely than the one 
to produce a single one, whatever is x. 

(iv) Whenever a is such that — G (a) has no constant sign, then both killing and 
branching can simultaneously occur at the death of the mother particle. A (x) may 
be put under the form A (x) — b (x) — d (x) where b (x) and d (x) are the birth 
(branching) and death (killing) components of A (x) . 

(v) Suppose A (x) is bounded below and let A* = ^infxefo.i] ^{x) > 0. Then one 
may view A as 

A (x) == A, (/i (x) - 1) , 

where /i (x) > 0. In this case, branching occurs at rate A*. When the mother 
particle dies, it gives birth to a spatially dependent random number M (x) of par- 
ticles (where M (x) — A(/i(x)) and A{p,{x)) is a geometrically distributed ran- 
dom variable on {0,1,2,...} with mean p.{x) = 1 + A(x)/A,). With Pm{x) — 
P (M (x) ^m) ^pq (x) q (x)", m > 0, po (x) == 1+^, Qo (x) = I - Po (x) 

A (x) = A, ^ mprn (x) - 1 = A, ^ (to - 1) p,n (x) - po (x) 

\m>l J \m>2 y 

Thus, the decomposition A (x) — b (x) — d (x) holds, where b and d can be read from 
A (x) = A, ' ' 



1 + ,u (x) 1 + ,u (x) 



{vi) In some other examples, the killing/branching rate X — G (a) /a is bounded 
above and below. Then A may be put under the form 

A(x) = A* {p(x) - 1) , 

where A* = sup^^grp j^i |A (x)| and < /i (x) < 2. In this case, we can assume that 
M (x) can only take the values or 2 with probability po {x) and p2 (x) respectively, 
with Po {x) +P2 (x) = 1. Then, /x (x) = E (M (x)) = 2p2 (x) and 

A (x) = A* (2p2 (x) - 1) = A* ip2 (x) - Po (x)) , 

giving a simple decomposition of A in the form A (x) = b (x) — d (x) with the mother 
particle living a random exponential time now with constant rate A, before giving 
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birth to none or two descending particles (a binary branching process). Note that 
Po (x) > P2 (x) (respectively po (x) > p2 {x)) when /i (x) < 1 (/i (x) > 1). 

Examples of a. When {xt;t > 0) is non-conservative, consider 

a(x)=Ef / c{xs)ds + d{xrjj , 

where the functions c and d are both assumed nonnegative on / and dl ~ {0, 1}. 
Then a > solves the Dirichlet equation —Ga{x) = c{x) > on / (= d{x) on 
dl) and so a is super-harmonic or excessive. We refer to [10) for examples of Doob 
transforms based on such super-harmonic functions allowing to understand various 
conditionings of interest when the starting point process {xt;t > 0) is a neutral WF 
diffusion or a WF diffusion with selection. 

Whenever a is super-harmonic for G, then /3 = l/a > is sub-harmonic for G = 
G + —g^dx- This results from the obvious identity 

showing that —Ga > entails —G (/3) < 0. 
Whenever [xt ]t > 0) is conservative and ergodic 

-Ej; / c{xs)ds -^ fJ-{c).= / c{y)d^i{y) 
t Jo *^°° Jo 

where dfi — m (y) dy is the invariant probability measure of {xt; t > 0) . Define 

ft 



lim Ex / c (xg) ds — tji (c) = a {x) . 



Thus 



a{x):^l {Ex{c{xs))- ^i,{c))ds 
Jo 

solves the Poisson equation 

—Ga (x) — c{x) :— c (x) ~ ^i{c) . 
We conclude that a is G— super-harmonic if ever c (x) > ^j,{c), \/x. o 

Background (multiplicative functional and path integral). The Doob trans- 
forms used here are particular instances of more general transformations based 
on multiplicative functionals. Let xt be the diffusion process ([ij governed by 
G = fdx + \9^dl with xq = x. 

Define the multiplicative functional Mt as the solution of the differential equation 

dMt = Mt ■ (a {xt) dt + 6 {xt) dwt) , Afo = 1, 
where a and b are arbitrary twice differentiable functions. Integrating, we get 
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Let _B be a Borel subset of /. Define a new process whose density p is obtained 
after a modification of the original one while using the multiplicative modulation 
factor Mt as 

/ p (x; t, y) dy := E^ [Mjl (x* e B)] = f E^ [Mt \xt^y]p (x; t, y) dy. 

J B J B 

Integrating Mt over paths with fixed two endpoints x and y, E^, [Mt [ Xt — y\ can 
be interpreted as the Radon- Nykodym derivative of p with respect to p, the density 
of Xt ■ By duality, let 

V (x, t) ^ E, [MtV (xt)] ,v{x,Q)^4, (x) . 
Applying Ito calculus, we get 

dtv = G (w) = (G + gbd^ + a) {v) =: (g + a) {v) , 

where the modified backward infinitesimal generator G is obtained by adding a 
drift term gbdx to G to produce G and a multiplicative part a-. The adjoint KFE 
giving the evolution of p is thus 

dtp ^ G* (p) = [G* + aj (p) , p (x; 0, y) = Sy (x) . 



- (Cameron-Martin-Girsanov) For instance, when a — and b = —f/g, the gener- 
ator of the transformed diffusion is G = \g^d^ killing the drift term of the original 
process governed by G. In this case, 



Mi 



^ -|/o[(f)'(-»)rf-+2f(-e)'i-=] 



Clearly in this case Mt is a martingale with 'Ex {Mt) — 1, assuming b to be bounded. 
This construction kills the drift of the original process while using a change of mea- 
sure. 

- (Feynman-Kac) When 5 = 0, the generator of the transformed diffusion is G = 
G + a adding a multiplicative component a to the one G governing the original 
process. In this case 

Mt = e^» a(x,)ds 

is the exponential of the integrated rate. liv{x,t) — E^ [Mttp (xt)], v (x,0) — tp (x) , 
then V solves 

dtv ^G{v)^{G + a) {v) , w (x, 0) = V (x) . 
In particular, if v (x, t) ~ Ex [Mt], v (x, 0) = 1 (x G (0, 1)) , then v solves 
dtv = G (w) = (G + a) {v) , u (x, 0) = 1 (x e (0, 1)) . 



(Doob) Suppose now 

dMt = Mt (a'l [xt) da (x*)) , Mq = 1. 
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This Mt is a particular instance of the general Mt introduced above. Indeed, 
applying Ito calculus, 

a~^ {xt) da {xt) — aT^Oi [fdt + gdw] + -a^^a"g^dt, 



leading to 



2 



a-^ ( fa' + ^a" ] ^ G {a) /a =: A (x) 



2 

b = q;^ a'g. 
Thus G = G + gbdx + a = G ^ a~^ a' g^ dx + A {x) as already observed earlier. 
Now, from the differential generation of Mj, 

a [x) 

only depends on the terminal and initial values of (xs ; < s < t) and not on its inter- 
mediate values (such a particular Doob transformation is thus a gauge). Thus here 

E^ [Mt \ Xt = y] = ^T^l consistently with the definition p (x; t, y) = £®p {x; t, y) . 

A super-harmonic example. 

Although this work chiefly focuses on Doob-transforms where branching is present 
in A, let us give a significant example where the Doob transform just produces 
killing like in (i). Suppose {xt]t > 0) is a non-conservative diffusion. Let Ai be 
the smallest non-null eigenvalue of the infinitesimal generator G of (xt ; t > 0) . Let 
a = Ml be the corresponding eigenvector, that is satisfying —Gui = XiUi > with 
boundary conditions ui (0) = ui (1) = 0. Then c = AiMi. The new KB operator 
associated to the transformed process {xt;t > 0) is 

(14) G(-) = -G(a).+G(-) = -Ai-+G(-), 

a 

obtained while killing the sample paths of the process {xt;t > 0) governed by G at 
constant death rate d = Xi. The transition probability of the transformed stochastic 
process {xt;t > 0) is 

P [x; t, y) = \ . p (x; t, y) . 
Ul (x) 

Define p{x;t,y) = e^^*p{x]t,y) . It is the transition probability of the process 
{xt ; t > 0) governed by G; it corresponds to the original process {xt ; i > 0) condi- 
tioned on never hitting the boundaries {0, 1} (the so-called Q— process of [xt] t > 0), 
see [H]). It is simply obtained from {xt;t > 0) by adding the additional drift term 

—g'^ to /, where Ui is the eigenvector of G associated to its smallest non-null 
eigenvalue. The determination of a = ui is a Sturm-Liouville problem. When t is 
large, to the dominant order 

Hx;t,y)^e-^*!^iM^, 

{ui,Vi) 

where vi is the Yaglom limit law of {xt;t > 0) . Therefore 

(15) pix;t,y) ^ ^^.t^^-x.t uiix)v,iy) ^ u,{y)v,{y) _ 

Ui{x) {Ui,Vi) {Ui,Vi) 



BRANCHING DIFFUSION AND SELECTION 15 

Thus the hmit law of the Q— process (xt; i > 0) is the normahzed Hadamard product 
of the eigenvectors ui and vi associated respectively to G and G* . On the other 
hand, the limit law of (xt', t > 0) is directly given by 

r^a\ -/ ^ \ ~l \ 1 2/« \,\ ' dz u{{y) 2fyJ^dz 

t^oc Zg'^ (y) Zg^ {y} 

where Z is the appropriate normalizing constant. Comparing (jlSp and ([1 



/ X ui{y) 2fyJ^dz , . , . 

vi{y) = ^-—re -"^ ^^(^ =ui{y)m{y). 

9 iy) 

The eigenvector vi associated to G* is therefore equal to the eigenvector ui associ- 
ated to G times the speed density of {xt;t > 0) . 

When dealing for example with the neutral WF diffusion (see Section 4 for ad- 
ditional details), it is known that Ai = 1 with ui — x (1 — x) and vi = 1. The 
Q— process {xt;t > 0) in this case obeys 



(17) dxt = (1 - Sxt) dt + ^Jxt (1 - xt)dwt, 

with an additional stabilizing drift toward 1/2: / (x) — —g^ (a;) = 1 — 2x. 

The limit law of the Q— process {xt;t > 0) in this case is 6j/(l — y). The latter 
conditioning is more stringent than the Yaglom conditioning and so the limiting 
law has more mass away from the boundaries (compare with the uniform quasi- 
stationary Yaglom limit ([8]) with vi = 1). 

3. The Wright-Fisher example 

In this Section, we briefly and informally recall that the celebrated WF diffusion 
process with or without a drift may be viewed as a scaling limit of a simple two 
alleles discrete space-time branching process preserving the total number N of in- 
dividuals in the subsequent generations (see |13j . [7], for example). 

3.1. The neutral Wright-Fisher model. Consider a discrete-time Galton Wat- 
son branching process preserving the total number of individuals in each gener- 
ation. We start with N individuals. The initial reproduction law is defined as 

follows: Let |kjv| := J2m=i^"^ ~ ^ ^^^ ^^ ■" (^ii-'-i^Af) be integers. Assume 
the first-generation random offspring numbers v>n := {vn {^) , ■■■,'^n {N)) admit 
the following joint exchangeable polynomial distribution on the discrete simplex 

|kjv| =N: 

This distribution can be obtained by conditioning N independent Poisson dis- 
tributed random variables on summing to A^. Assume subsequent iterations of 
this reproduction law are independent so that the population is with constant size 
for all generations. 

Let Nr (n) be the offspring number of the n first individuals at the discrete gener- 
ation r e No corresponding to (say) allele Ai (the remaining number N — Nr (n) 
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counts the number of alleles A2 at generation r). This sibship process is a discrete- 
time Markov chain with binomial transition probability given by: 



P (iV.+i (n) = k'\Nr (n) = k)= {^)j {j^^ 



1 _ h: "''' 



Assume next that n — [Nx] where x £ (0, 1) . Then, as well-known, the dynamics 
of the continuous space-time re-scaled process Xt '■— N[pft] ('^) /N, t G R_|- can 
be approximated for large N , to the leading term in N^^ , by a Wright-Fisher-Ito 
diffusion on [0, 1] (the purely random genetic drift case): 



(18) dxt = \/xt (1 - xt)dwt, xq = X. 

Here {wt\t > 0) is a standard Wiener process. For this scaling limit process, a unit 
laps of time t = 1 corresponds to a laps of time N for the original discrete-time 
process; thus time is measured in units of N . If the initial condition is a; = N~^ , xt 
is the diffusion approximation of the offspring frequency of a singleton at generation 

[Nt]. 

Equation (fT8| is a 1— dimensional diffusion as in ([T]) on [0, 1] , with zero drift / (x) = 
and volatility g (x) = -^/a; (1 — x). This diffusion is already in natural coordinate 
and so if (x) = x. The scale function is x and the speed measure [x (1 — x)] dx. 
One can check that both boundaries are exit in this case: The stopping time is 
Tx — Txfi A Tx^i where t^^q is the extinction time and t^.i the fixation time. The 
corresponding infinitesimal generators are G {■) — ^x{l — x)dl{-) and G* {■) = 

3.2. Non-neutral cases. Two alleles WF models (with non-null drifts) are clas- 
sically obtained by considering the binomial transition probabilities hm{N,pN) : 

P (Nr+i {n) ^k'\Nr in) ^ k) = {^^)j {pN {^^y^ {l-pN (^^' ' 

where 

pjv (a;) : x e (0, 1) ^ (0, 1) 

is now some state-dependent probability (which is different from the identity a;) 
reflecting some deterministic evolutionary drift from the allele Ai to the allele A2. 
For each r, we have 

E {Nr+i (n) I Nr (n) = k) = NpN Q 
a^ {Nr+i (n) I Nr {n) ^ k) = Np^ (^^ (^1 - Pn (J^ 

which is amenable to a diffusion approximation in terms of Xt '■— -/VfArti (n) /N, 
t € R+ under suitable conditions. 

For instance, taking 

Pn (x) = (1 - 772, at) X + TTl^N (1 - x) 
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where (7ri^Ar,7r2,Ar) are small (A'^- dependent) mutation probabilities from A2 to Ai 
(respectively Ai to A2). Assuming {N ■ tti^n,N ■ tt2,n) -^ (7ri,7r2), leads after 

scaling to the drift of WF model with positive mutations rates (tti, 7r2). 
Taking 

(1 + si,n)x 

PN [X) = — — r 

1 + Si,nX + S2,N (1 - X) 

where s; at > are small A^— dependent selection parameter satisfying N -Si n ^ 

(Ti > 0, i — 1, 2, leads, after scaling, to the WF model with selective drift / (x) = 
ax {1 — x), where a :— ai—a2- Typically, the drift / (x) is a large N approximation 
of the bias: N {p^ (x) — x) . The WF diffusion with selection is thus: 



(19) dxt = axt (1 - xt) dt + ^J xt (1 - Xt)dwt 

where time is measured in units of N. Letting 9t — Nt define a new time-scale with 
inverse te — 6/N, the time-changed process ye ~ xgjf^ now obeys the SDK 



dye = sye (1 - ye) dO + \l j^Ve (1 - ye)dwe, 

with a small diffusion term. Here s = si ^ S2 and time 9 is the usual time-clock. 

The WF diffusion with selection ^^ tends to drift to o = 1 (respectively o = 0) 
if allele Ai is selectively advantageous over A2 : ai > 02 (respectively a\ < 02) in 
the following sense: if cr > (respectively < 0), the fixation probability at o = 1, 
which is [16] 

1 — e^^*^^ 
P(r,,i<^-,o) = Y3— ^, 

increases (decreases) with a taking larger (smaller) values. 

The usual way to look at the WF diffusion with mutation and selection is to com- 
pose the two above mechanisms ptq (x) corresponding to mutation and selection 
respectively. In the scaling limit, one obtains the standard WF diffusion model 
including mutations and selection as: 



(20) dxt = [(tti - (tti + 7r2) xt) -f axt (1 - xt)\ dt + \J xt (1 - Xt)dwt. 



4. The neutral WF model 

In this Section, we particularize the general ideas developed in the introductory Sec- 
tion 2 to the neutral WF diffusion (|18l) and draw some straightforward conclusions 
most of which are known which illustrate the use of Doob transforms. 

4.1. Explicit solutions of the neutral KFE. As shown by Kimura in (fT5|), the 

Kolmogorov forward (and backward) equation is exactly solvable in this case, using 
spectral theory. The solutions involve a series expansion in terms of eigen-functions 
of the KB infinitesimal generator with discrete eigenvalues spectrum. 

Let A/c = k (k + 1) /2, k > 0. There exist Uk = Uk{x) and Vk — Vk{y) solving 
the eigenvalue problem: —G{uk) = XkUk and —G* (vt) = XkVk- With {vk,Uk) — 
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/o Uk {x) Vk {x) dx, the transition probability density p {x; t, y) of the neutral WF 
models admits the spectral expansion 

p {x; t,y) = y^ bke~^''*Uk (x) Vk (y) where bk = -. r 

where Uk (x) are the Gegenbauer polynomials rescaled on [0, 1] and normalized 
to have value 1 at a; = 0. In particular, uq [x) = x, ui (x) = x — x"^, U2 {x) = 
x-ix'^ + 2x^, Us, (x) = x-6x'^ + l0x^-5x^, m (x) = x-10x'^+30x^-35x'^ + Ux^, ... 

Next, Vk (y) = m {y) Uk (y) where m {y) = 1/ [y (1 — y)) is the speed density of the 
neutral WF diffusion. For instance, wo (y) — j-^i '^i iv) — 1: ^2 (y) = 1 — 2y, 
^3 (y) = 1 - 5y + 5y2, V4, (y) = 1 - 9y + 21y^ - Uy^,... 

Although Ao = really constitutes an eigenvalue, only vq (y) is not a polynomial 
and the spectral expansion of p should start at k = 1, expressing that p is a sub- 
probability. When fc > 1, from their definition, the Uk (x) polynomials satisfy 
Uk (0) = Uk (1) = in such a way that Vk (y) = m (y) • Uk (y) , k > 1 is a polynomial 
with degree fc — 1. 

The series expansion for p {x; t, y) solves the KFE of the WF model. 
We have P (tj. > f ) = J„ P {xf G dy) and so 

p, ix) := P (r. > t) = 5: &P^e-^^W (x) 

is the exact tail distribution of the absorption time. 

Since vi (y) = 1, to the leading order in t, for large time 

P {xt e dy) = 6e"* ■ x {I ~ x) dy + O (e^^*) 

which is independent of y. Integrating over y, p^ (x) :— P (t^ > t) ^ 6e^* -x (1 — x) 
so that the conditional probability 

(21) P {xf edy\T,>t) - dy 

t^oo 

is asymptotically uniform in the Yaglom limit. As time passes by, given absorption 

did not occur in the past, Xt — ^ Xoo (as t — J> 00) which is a uniformly distributed 
random variable on [0, 1]. 

4.2. Additive functionals for the neutral WF and Doob transforms. Let 

{xt\t > 0) be the neutral WF diffusion model defined by (fT8|) on the interval / = 
[0, 1] where both endpoints are absorbing (exit). Consider the additive quantities 



a{x)='E( c{xs)ds + d{xrj] 



where functions c and d are both nonnegative. With G = ^x (l — x) d"^, a{x) 
solves: 

-G{a) = ciixel 

a = d if X e dl. 

Therefore a is a super-harmonic function for G. 
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Take c = lini£4.o ^l(j,_e_y+£) (x) =: 6y (x) and d ~ 0, when y G / : in this case, 
a :— Q{x,y) is the Green function (the mean local time at y given the process 
started at x). The solution takes the simple form 

X 

3ix,y) = 2- iix <y 

y 

I ~ X 
3ix,y) = 2- if a; > y. 

The Green function solves the above general problem of evaluating additive func- 
tionals a (x) : 



a{x) = Q{x,y)c{y)dy iix e I 
a ~ d ii X £ dl. 

There are many interesting choices of c therefore leading to a, allowing to compute 
for example the mean time till absorption for the neutral WF diffusion, the prob- 
ability to hit state 1 before 0... For each choice of a, it is interesting to study the 
transformed process {xt', t > 0) whose transition probability is given by 

P{x]t,y) = —--p{x;t,y), 
a (x) 

in terms of the original process transition probability p (x; t, y) . This allows for ex- 
ample to understand the neutral WF process conditioned on exit at some boundary 
and to evaluate for this new process interesting average additive functionals such 
as the mean time needed to hit the exit boundary... For detailed similar examples 
arising in the context of WF diffusions and related ones, see [TU] . 



5. The WF model with selection 

Now we briefly focus on the diffusion process P^ . Let {vu {y))k>i ^^ the Gegen- 
bauer eigen-polynomials of the KF operator corresponding to the neutral WF diffu- 
sion (US]), so with eigenvalues Afc = fe (fc -I- 1) /2, fc > 1. Define the oblate spheroidal 
wave functions on [0,1] as 

(22) <(y)-E/^'^'(y)' 

i>i 

where /^ obey the three-term recurrence defined in [5D]. In the latter equality, the 
I summation is over odd (even) values if k is even (odd) . 

Define v% (y) = e'^^w^ (y) and u^ (x) = i^^j^v"^ (x) where m (x) = e^"'-^ j {x (1 - x)) 
is the speed measure density of the WF model with selection (fTO]). 

The system (u^ (x) ,w^ {xS)h>\ constitute a system of eigen-functions for the WF 
with selection generators —G and —G* with eigenvalues A^ implicitly defined in [20) . 
thus with —G (u'j!,) = A^u^ and —G* (w^) = A^w^. The eigen-function expansion of 
the transition probability density of the WF model with selection is thus, |14) : 



(23) p (x; i, 2/) = V ble-^'-'ul (x) v^, (y) 



k>l 
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where 6^ — (f^,""^) . The WF model with selection can be viewed as a per- 
turbation problem of the neutral WF model (see [Sj). There exist perturbation 
developments of A^ around X^ with respect to a'^, [14! . They are valid and useful 
for small a. 



The WF diffusion process Xt with selection (|T9|) is non-conservative, with finite 
hitting time Tx of one of the boundaries. Following the general arguments developed 
in Section 2, the Yaglom limit of Xt conditioned on r^; > i is the normalized version 
of 

(24) < (y) = e'^yw-, (y) . 

The limit law of xt conditioned on never hitting the boundaries in the remote future 
is the normalized version of 

(25) < (y) «^ (y) = ^< {xf = y (1 - yX {yf . 

m(y) 

Because the latter conditioning is more stringent than the former, the probability 
mass of ((25|) is more concentrated inside the interval than ([24]) . Compare with the 
statements at the end of Section 2 concerning the neutral WF diffusion. 



6. DOOB TRANSFORM OF THE NEUTRAL WF MODEL: SUB-CRITICAL BD 

In this Section, we define the branching WF diffusion model with selection while 
applying a Doob transform to the neutral WF model, based on the sub-harmonic 
additive functional a (x) = e°'^, say with a > 0. We then study in detail the ob- 
tained branching process. 



The starting point is thus the neutral WF diffusion: dxt — \/ Xt (1 — xt)dwt, xq = 
x€ (0,1). 

For this model, G = ^x {I — x) d^ and both boundaries are exit. With Xf. = 
k {k + 1) /2, fc > 0, its transition density p {x; t, y) admits the spectral representa- 
tion 

(26) p {x; t,y)^Yl ^ke-^^'uk {x) Vk {y) , 

fc>i 

in terms of the Gegenbauer eigen-polynomials (see Subsection 4.1). We shall 
consider the following transformation of paths on the neutral WF model: Let 
a{x) = e*^^, cr > and consider G {■) = a^^G{a-) — G{-) + b{x) ■ . We now 
have G (a) = \a'^x (1 - x) e"^^ and so h{x)^G (a) ja ^^x(\~x)> 0. 
Note that —G (a) < indicating that a is sub-harmonic for G. 

In this case study, one selects sample paths of {xt ; i > 0) with large a (y) and we 
claim that this is an alternative interesting way to introduce selection in the neutral 
WF diffusion process. 

The dynamics of {xt\t > 0) governed by G is easily seen to be the standard WF 
with selection dynamics ([19]) 



dxt = axt (1 - Xt) dt + y/xt (1 - Xt)d' 



Wt, 
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subject to additional quadratic branching at rate b (x) = \cr'^x (1 — x) inside /. We 
indeed have 



G(-) = e-'^"G(e""-) = fc(a;)-+G(-), 
where 

G =: /a, + hfdl - ax (1 - x) d, + \x{l~ x) dl 

is the KBE operator of the dynamics (xt ; i > 0) . Recall that xt is transient and so 
hits one of the boundaries {0, 1} in finite time Tx- 

To summarize, in our branching diffusion way to look at selection, we move from 
the neutral WF diffusion {xt]t > 0) to the standard WF diffusion with selection 
{xt;t > 0) but subject to additional branching at rate b {x) . 

Remark. With (3 (x) :— a (x) = e^"'^ , we clearly have 

G(/3(x)) = 

and /3 is an harmonic function for G and as a result, Doob-transforming G by /3, 
we get 

which is the infinitesimal generator of the original neutral WF martingale, o 



/3-1g(/3.) = M)^'G(«/3-)-G(-) 



The birth (creating) rate & > in G is bounded from above on (0, 1). It may be 

2 

put into the canonical form b (x) — 6* (u (x) — 1) where 6* = max ib (a;)) = ^ > 

2:e[0,l] 

and 

(27) m(-t) = l+4a;(I-a;), 

whose range is the interval [1, 2] as a; £ [0, 1] . 

The density of the transformed process is p{x;t,y) = ^t|tp (x; t, y) . It is exactly 
known because so is p is from (|26p . 

The transformed process (with infinitesimal backward generator G) accounts for 
a branching diffusion (BD) where a diffusing mother particle (with generator G 
and started at x) lives a random exponential time with constant rate 6* . When the 
mother particle dies, it gives birth to a spatially dependent random number M (x) 
of particles (with mean /i(a;)). M (x) independent daughter particles are started 
where their mother particle died; they move along a WF diffusion with selection 
and reproduce, independently and so on. 

Because /i (x) is bounded above by 2 and larger than 1 (indicating a super-critical 
branching process), we actually get a BD with binary scission whose random off- 
spring number satisfies ('w.p.' meaning 'with probability') 

M (x) = w.p. po (x) = 

M{x) = 1 w.p. pi{x) ^2- fi (x) 

M (x) = 2 w.p. p2 (x) = /i (x) — 1, 

with p2 (x) > pi (x) (the event that 2 particles are generated in a splitting event is 
more probable than a single one). 
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For such a transformed process, the trade-off is as follows: there is a competition 
between the boundaries {0, 1} which are absorbing for the particle system and the 
number of particles Nt (x) in the system at each time t, which may grow due to 
binary branching events (or remain steady when M (x) = 1). 

The density p of the transformed process has the following interpretation 



(28) 



p{x;t,y) = E 



Nt{x) 



where p'"-' (a;; i, y) is the density at (i, y) of the nth alive particle descending from 
the ancestral one (Eve), started at x. In the latter formula, the sum vanishes if 
Nt {x) = 0. A particle is alive at time t if it came to birth before t and has not been 
yet absorbed by the boundaries. 

Let P( {x) — J/p j^N p (x; t, y) dy. Then pj (x) is the expected number of particle alive 
at time t. We have 

dt-pt (x) ^Gipt (x)) , po (x) ^l{xe (0, 1)) . 

Remark. From the Feynman-Kac formula, p in (l^5|) is also 

p {x- 1, y) = E, (e^"^'"= "(^^^''^ \xt=y)p {x; t, y) 

and 

p,(x)=E.(er'^''(^=)'^^).o 

But then g (x; t, y) :— p {x; t, y) /j)f. (x) obeys the forward PDE 

dtPt (x) 



dtq {x; t,y) = - 



Ptix) 



+ biy)]q{x;t,y) + G*{qix;t,y)) 



as a result of dtp (x; t,y) — G [p (x; i, y)). We have 



(29) 



qix;t,y) 



E 



En'J^^P^'^Hx-Ay) 



E[Ntix)] 

showing that q (x; t, y) is the average presence density at (i, y) of the system of 
particles all descending from Eve started at x. 



Clearly - 
because 



log Pt(^) 



-> Ai 



1 (and therefore also — -^(x) ^^ ^^ Hospital rule), 



Pt{x) 



—r-r^bke ^"'ufe (x) / a{y)vk{y)dy. 
a (x) ^—' 



k>l 



The expected number of particles in the system decays globally and exponentially 
at rate Ai. 

The BD transformed process therefore admits an integrable Yaglom limit q^, so- 
lution to -G* [q^) = (Ai +b{y))q^ or ~G* [q^) = Xiq^. With vi (y) = 1, the 
first eigenvector of — G* associated to the smallest positive eigenvalue Ai = 1, g^^ 
is of the product form 

(30) q^{y)^C.e-yv,{y) = -^. 
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The arbitrary multiplicative constant C* was chosen in such a way that q^ (y) is a 
probability. 

By analogy with the Yaglom construction, this limiting probability q^ca.n be called 
the quasi-stationary Yaglom average density at (i, y) for the BD particle system (it 
is also the ground state for G ). 

There is also a natural eigenvector cf)^ of the backward operator — G, satisfying 
— G {(j)oo) — -^1000 (the ground state for G). It is expHcitly here 

(31) ^^(:c)^-^u,{x) = ^^^—^e-'^^x{l-x). 

a (x) a 

The arbitrary multiplicative constant G = 6/G* was chosen in such a way that 
/o 9oo (y) 4>oo (y) dy — I. Note that the spectral structures of both G and G are 
easily obtainable from the ones of G* and G thanks to the Doob transform structure. 

In the terminology of [22], both operators G (•) + Ai- and its adjoint are critical O 
In this context, the constant Ai is called the generalized principal eigenvalue. The 
eigen-functions ((/'ooi 9oo) ^^'^ their associated ground states. We note that we have 
the L^— product property (See [22], Subsection 4.9). 

i_ _ /-i 

4'oo (^) ^oo (^) dx = 6 ui {x) vi {x) dx = \ < oo. 
Jo 



Remark. Using the Fcynman-Kac representation of pj (a;) , we get 
-- logE^ fe-^'o"^" b{x,)ds\ _^ ^^ ^ 1 and 



With p„ (x) = P (M (x) = m), let 

l{x) = ^^ Prn (x) m log m = 2 log 2p2 (x) . 

m>l 

We have the x log x condition: 

(32) / l{x) (j)^{x)q^{x)dx — 48\og2 x {1 — x)ui{x) vi{x) dx < oo. 

Jo Jo 

We conclude (following jlj and [2J) that, as a result of the condition ([5^ being 

trivially satisfied, global extinction holds in the following sense: 

(i) P (Nt (x) = 0) -^ 1, uniformly in x. 

(a) there exists a constant 7 > : e'^i* [1 — P {Nt (x) = 0)] — > 7</>oo (2;) , uniformly 
in X. 



G()-|-Ai- (G (■)+Ai-) is said to be critical if there exists some function (/);jQ G C^ (respectively 
9oc G C^), strictly positive in (0,1) , such that: G (<t>ao) + ^I'f'oo — (respectively G {q^ + 
Ai5;,„ = 0) and the operators do not possess a minimal positive Green function. 
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(iii) For all bounded measurable function ij) on I : 

'Nt(x) 

^^(5?("))|iV,(x)>0 



E 



*^°° 7(0,1) 



From (i), it is clear that the process gets ultimately extinct with probability 1. In 
the trade-off between pure branching and absorption at the boundaries, all particles 
get eventually absorbed and the global BD process turns out be sub-critical (even 
though /i (x) = EM (x) > 1 for all x G (0, 1)): Probability mass escapes out of / 
although the BD survives with positive probability. 

In the statement (m) , the quantity 1 — P {Nt [x) = 0) = P {Nt (x) > 0) is also 
P (r {x) > t) where T {x) is the global extinction time of the particle system de- 
scending from an Eve particle started at x. The number — Ai is the usual Malthus 
exponential decay rate parameter. From {ii) , (j)^ (x) has a natural interpretation 
in terms of the propensity of the particle system to survive to its extinction fate: 
the so-called reproductive value in demography. 

[iii) with -0 = 1 reads E [Nt (x) | Nt (x) > 0] — ?► ^~^ giving an interpretation of 

t—^oo 

the constant 7 (which may be hard to evaluate in practise). 

The ground states of G -t- Ai and its adjoint are thus (0oo;9oo) ^'^'^ explicit here. 
It is useful to consider the process whose infinitesimal generator is given by the 
Doob-transform 

TJ (G -f Ai) (0^-) = 0~' (g + 6 + Ai) (0^-) , 

because product-criticality is preserved under this transformation. The ground 
states associated to this new operator and its dual are (l, (?!>oo9oo)- Developing, we 
obtain a process whose infinitesimal generator is 

with no multiplicative part. In our case study, we get ■^x (1 — x) 9^ + (1 — 2x) dx 
adding a stabilizing drift towards 1/2 to the original neutral WF model. The asso- 
ciated diffusion process is positive recurrent and so its invariant measure 0oo9oo = 
6uiVi — 62/(1 — y) is integrable with mass 1. It is the beta(2,2) limit law of the 
Q— process (see (fT7|) and the comments at the end of Section 2 relative to the neu- 
tral WF diffusion). 

Remarks. 

/, ,\Nt{x) 

(i) At time t, let I x^ j denote the positions of the BD particle system. Let 

V / n—l 



u (x, t] z) — 'El 



n 



Wt(:r)^ (£{"') 



stand for the functional generating function (|2:| < 

1) of the measure-valued branching particle system, u (x, t; z) obeys the nonlinear 
(quadratic) Kolmogorov-Petrovsky-Piscounoff PDF, [17j : 

dtu (x, t; z) = b^9 (x, u (x, t; z)) + G {u (x, t; z)) ; u (x, 0; z) = z^^^' , 
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where 9 (x, z) = E [z^^^^'] — z = (p2 (x) z^ + pi (x) z) — z or 

(x, z) = 4x (1 — x) z (z — 1) 

is the shifted probabihty generating function of the branching law of M (x) . Thus, 
the nonhnear part reads b<,6 (x, u (x, t; z)) — h (x) u (x, i; z) (u (x, t; z) — 1) which is 
quadratic in u. 

In particular, if u (x, t) := dzU (x, t; z)^^^ — E X^nli V" ( ^t ) i "" (^^i obeys the 
linear backward PDE 

dtu (x, t) — b (x) u (x, t) + G {u (x, t)) ; u (x, 0) = i/i (x) 
involving G (■) = G (•) + b (x) •. We have the Feynman-Kac interpretation 

The latter evolution equation is the backward version of the forward PDE giving 
the evolution of p (x; i, y) as dtP (x; t,y) = G {p (x; t,y)) ,p (x; 0, y) = S^ (y) . 

(m) Let us look at the branching diffusion process governed by G would time be 
measured using the time substitution 

t nt 

Jo 

for each of the particles that came to birth before t. 

Then G ^ Q := ^(i_^) G = adx + ■^d^ + ^cr^ ■ . In particular, each motion yg = Xtg 

is a Brownian motion with constant drift (a Gaussian process). This new Q is the 
one of absorbing Brownian motion with drift a on [0,1], including branching at 
constant rate hcr^ ■ The Sturm-Liouville problem for Q admits the eigenvalues A^ = 



2 



, A: > 1 with eigen-states Uk (x) (x e "^ sin [kirx) and Vk (y) oc e'''^ sin (kiry) . 

The spectral gap is Ai — ^ > ^'^'^ ^^^ time-changed branching diffusion 
also becomes eventually extinct, sub-critically: The time substitution changes the 
spectral structure of the model but not its qualitative features, o 



7. DOOB TRANSFORM OF THE WF MODEL WITH MUTATIONS: CRITICAL BD 

In this Section, we start from the WF model with mutations. Using the same Doob 
transform based on the additive functional a (x) = e'^^ to introduce selection, we 
end up with a WF diffusion process with killing and branching describing the effect 
of selection on the WF model in the presence of mutations. We show that in this 
setup, the resulting branching diffusion process is no longer sub-critical; rather, it 
turns out to be critical. 

Suppose the starting point model is now the WF diffusion with mutations: 



dxt — (tti — TTXt) dt + \J Xt (1 — Xt)dwt,xo = x G (0, 1) , 
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with TT := TTi + ■K2- For this model, G = (tti — nx)dx + -^x {1 ~ x) d^ and both 



- 'Kx)d 
tingfl. 



value 1 at X = 0. In particular, mq (x) = 1, ui (x) = 1 — —x, U2 (x) = 1 — —x 



boundaries are chosen as being entrance (reflectingjj. The WF diffusion process 
with mutations is now ergodic. With 

kjk-l+TT) 
^k = ^ , fc > 0, 

its transition density p (x; t, y) now admits the discrete spectral representation 
(33) p {x; t,y) ^Yl bke^^^^k {x) Vk (y) ■ 

fe>0 

Here, Uk (x) are the Jacobi polynomials rescaled on [0, 1] and normalized to have 
value 1 at X = 0. In particular, uq (x) = 1, ui (x) = 

^tta+^t)^ ^''- ^«^t' ^fc (y^ = "^ (y^ ""^ (y^ ^^^^-^ 

^^y' r(27ri)r(2^2)^ ^ ■'^' 

is the speed density of the ergodic WF diffusion with mutations (its normalized 
invariant measure). Note that the fc = term in ([33]) is precisely m (y) as required. 
Because the transition probability density of the WF diffusion with mutations has 
also a discrete spectral representation, this model is amenable to a similar analysis 
than the neutral WF diffusion. 



Proceeding as for the neutral case, we shall consider the following transformation 
of paths for the WF model with mutations: Let a (x) — e""^ and consider a trans- 
formed process with infinitesimal generator G {■) — a^^G{a-) . The multiplicative 
part of G is now 

2 

X{x) — G {a) /a — (J (tti — nx) + —x (1 — a;) . 

Note that now a is neither sub-harmonic nor super-harmonic for the infinitesimal 
generator G including mutations because the sign of A (x) varies as x varies. 

In this case study, one selects sample paths of the WF diffusion model with muta- 
tions {xt', t > 0) with large terminal values of a (y) . The dynamics of (xt] i > 0) is 
easily seen to be the WF with mutation and selection dynamics of the type (I^UI 

dxt = [(tti - TTXt) + axt (1 - xt)] dt + \/xt (1 - xt)dwt, 

subject to additional quadratic killing and branching at rate A (x) inside /. We 
indeed have 



G (•) = e-'^^'G (e^^-) = A (x) • +G (•) , 

where G = [(tti — nx) + crx (1 — x)] dx + ^x (1 — x) 9^ is the KBE operator of the 
dynamics (xt;i > 0) . 



When both the mutation rates mi and U2 are greater than 1/2, the boundaries are entrance. 
When either u\ or U2 is smaller than 1/2 the corresponding boundary is regular and one needs 
to specify whether it is reflecting or absorbing or a mixture of the two. We force here the regular 
boundaries to be entrance. 
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To summarize, in our branching diffusion way to look at the action of selection, we 
move from the WF diffusion with mutations {xt;t > 0) to the standard WF diffusion 
with mutation and selection [xt] t > 0) , but subject to additional killing/branching 
at rate A (x) . 

Remark. With /3 = «-! = e"'"^, again G(^) = and ^"^G(^-) = G {■) is the 
infinitesimal generator of the original WF model, now with mutations, o 

The birth (creating) and death (annihilating) rate A in G is bounded from above and 
below on (0, 1). It may now be put into the canonical form A {x) — A» {fi (x) — 1) 

where A, = max (|A(a;)|)and 

xe[o,i] 

(34) /,(:,) = i + ^ 

whose range belongs to the interval [0, 2] as a; £ [0, 1] . 

Note that when n > a/2, Xt, — a (tti V 112) whereas when n < a/2, A* = A (x*)V(77r2 
where x„ = 1/2 — ir/a > 0. 

The density of the transformed process is p{x;t,y) — ^r^p (x; t, y) . It is exactly 
known because p is known from (|33p . 

The transformed process (with infinitesimal backward generator G) accounts for a 
branching diffusion (BD) where a diffusing mother particle (with generator G and 
started at x) lives a random exponential time with constant rate A*. When the 
mother particle dies, it gives birth to a spatially dependent random number M (x) 
of particles (with mean /i (x)). If M (x) ^ 0, M (x) independent daughter particles 
are started where their mother particle died; they move along a WF diffusion with 
mutation and selection (with generator G) and reproduce independently, and so on. 

Because fj, (x) is bounded above by 2 and larger than 0, we actually get a BD with 
binary scission whose random offspring number satisfies 

M{x) =0 w.p. po (x) ^1- fi{x) /2 

M (x) = 1 w.p. pi (x) = 

M (x) = 2 w.p. p2 [x] = n (x) /2. 

Note that 

A (x) = A* (p2 (x) - Po [x)) =: b{x) - d (x) 
identifying the birth and death components of the full multiplicative rate A (x) . 

For such a transformed process, the trade-off is of a different nature: there is a 
competition between the boundaries {0, 1} which are now reflecting for the system 
of particles and the number of particles Nt (x) in the system at each time t, which 
may grow or diminish due either to branching or killing events. In the presence of 
mutations, the particles are no longer killed once they hit the boundaries, suggest- 
ing that there should be a greater amount of them alive in the system. However, 
in this new model, there is an opportunity to kill the particles inside the definition 
domain, when they branch. The question now being: does the new trade-off result 
in global extinction or global explosion of the particle system? We will now show 
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that critical global extinction occurs. 

The density p of the transformed process again has the interpretation (l28l) , where 
p^"^ {x;t,y) is the density at (t,y) of the nth alive particle descending from the 
ancestral one (Eve), started at x. In the latter formula, the sum vanishes if Nt (x) = 
0. A particle is alive at time t if it came to birth before t and has not yet been 
killed by a killing event. 

Let P( (x) = r.(^ -^^ p (x; t, y) dy. Then pj {x) is the expected number of particle alive 
at time t. We have 

dt-pt [x] ^G{pt (x)) , po (x) = 1 (x e (0, 1)) . 

But then q {x; t, y) :— p {x; t, y) j'p^ (x) obeys the forward PDE 

d,q{x-Xy)={~^-^4^+h{y)\q{x-t,y) + &{q[x-Xy)) 

as a result of 9ip {x\ t,y) — G {p {x; t, y)). We again have (1^ . with q {x; t, y) the 
average presence density at (t, y) of the system of particles all descending from Eve 
started at x. 

Clearly _-2I-£ti£i _> Aq = (and therefore also — L^,,' ), because 

Pt (x) = —TT y2 bke^^^^Uk (x) / a (y) Vk (y) dy. 
a (x) ^-^ Jr. 

^ ' A:>0 •'" 

The expected number of particles in the system decays globally at rate Ai towards 
the non-zero limiting value 

Poo (x) -"^ oi i^y^ bouo {x) a (y) vo {y) dy ^ e'""" e'^'" m [y) dy . 

Jo Jo 

The BD transformed process therefore admits an integrable Yaglom limit q^. so- 
lution to -G* (q^) ~ A {y)qaa or -G (q^o) = 0. With vo (y) = rn{y) , the first 
eigenvector of — G* associated to the smallest positive eigenvalue Aq = (the equi- 
librium density of the WF diffusion with mutations) , q^o i^ of the product form 

(35) goo (y) = 7T -— — ■ 

Jp CVm (y) dy 

This explicit limiting probability q^is the Yaglom limiting average presence density 
at {t,y) for the BD system of particles (it is also the ground state for G ). 

There is also a natural eigenvector cf)^ of the backward operator —G, satisfying 
— G {(jjoo) — (the ground state for G). It is exphcitly here 

(36) 0OO i^) = — TT^o (^) / ^"^^ (y) '^y = ^ '^'^ / ^'^^'^ (y) ^y- 

a [x) Jo Jo 

Both operators G (•) and its adjoint are again critical. The constant Aq = is the 
new generalized principal eigenvalue; The eigen-functions (^ooj^oo) ^■'"S the new 
associated ground states. 
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We note that we have the L^— product property 

/ uo (x) vo (x) dx = 4>^ {x) q^ {x) dx ^1 < oo. 
Jo Jo 

Clearly the ground states of — G and —G are defined up to arbitrary multiplicative 
constants. Note that we chose these constants in such a way that /„ q^ (y) dy = 1 
and /o ^oo (x) q^ (x) dx = l. 
With p„i (x) = P (M (x) = m), let 

K (x) = ^ m (m - l)pm (x) = 2p2 (x) . 

m>2 

Because p2 {x) is a degree two polynomial in x, we have the condition: 



K [x) (f)^ {x) q^ [x] dx — 2 p2 (x) uq (x) vq (x) dx < oo. 



(37) 



We conclude (following [Ij and |2J) that, as a result of the condition ([57]) being 
trivially satisfied, global extinction holds critically, in the following sense: 



(j) P {Nt (x) = 0) —T' 1, uniformly in x. 

{^^) Let /i* = E,t^i^ '^^j") , with f,, (^) = E^4^i^ ^ (^^"^) 
There exists a finite positive constant : 



"=sy. ■=■ 



Mt {<l>)^ ~ Mt (0^) ^oo {^) dx 



2i '' 



Mt (0)^ - A^t (0^) 



such that t [1 — P {Nt (x) = 0)] — > /i "'^^g^ (x) , uniformly in x. 
(mi) For all bounded measurable function ip on I : 

'Nt{x) 



t 



J2^[x[-^)\Nt{x)> 



,^ M / ^{y)qoo{y)dy. 

*^°° 7(0,1) 



From (i), it is clear that the process gets ultimately extinct with probability 1. In 
the trade-off between killing-branching and reflection at the boundaries, all particles 
get eventually absorbed but the global BD process turns out be critical. Thus, 
the killing part of A (x) is strong enough to avoid the explosion of the number of 
particles inside the unit interval, resulting in an overall critical process where global 
extinction still holds. 

In the statement (ii) , 1 - P (A^* (x) = 0) = P (T (x) > t) where T (x) is the global 
extinction time of the particle system. The Pareto tails of T (x) decay like t~^, thus 
algebraically slowly: the time till extinction in this critical model is much longer 
than in the previous neutral sub-critical case (with exponential tails). From (ii) , 
0OO (x) has again a natural interpretation in terms of the propensity of the particle 
system to survive to its extinction fate. 
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(Hi) with ip = 1 reads 7E [Nt (x) \ Nt {x) > 0] — > /i giving an interpretation of 
the constant fi. The constant fi is also ([2], page 287) 

H^-X., K{x)'4>^{xfq^{x)dx^Kj p2{.x)lj)^(xfq^{x)dx<GC, 
and so is exphcitly available in our case. 

The ground states of G + Aq and its adjoint are thus {(f)^ , q^) and explicit here. It 
is also useful to consider the process whose infinitesimal generator is given by the 
Doob-transform 

TJg i^^-) =TJ {g + \) (^^■) , 

because product-criticality is preserved under this transformation. The ground 
states associated to this new operator and its dual are (l, (j^aJloc)- Developing, we 
obtain a process whose infinitesimal generator is 

with no multiplicative part. The associated diffusion process is the starting point 
WF diffusion with mutations, which is positive recurrent and so its invariant mea- 
sure 0oo9oo = "^^o^o = fn (y) is integrable. 



Remark. The functional generating function u (x, t; z) = E 



|-rWt(x) ^i,[xl 



ln=l 



of 



the measure- valued branching particle system obeys now the nonlinear (quadratic) 
PDE: 

dtu{x,t;z) = \^,6{x,u{x,t;z)) + G{u{x,t;z)) ; u(a;,0;z) = z'^''^^ 

where 9 {x, z) — 'Ei [z*'-'*^^)] — z — (p2 {x) z^ + pa (a;)) — z or 

9 {x, z) = {z - 1) {p2 (x) z ~po (x)) 

is the shifted probability generating function of the branching law of M (x) . 

If u (x, t) := dzU (x, i; z)^^^ = E J2n=i^ V' f ^t" j , recalling A (x) = A* (p2 (x) - po (x)) , 
u {x, t) obeys the linear backward PDE 

dtu (a;, t) = \ {x) u {x, t) + G {u {x, t)) ; u{x,0) — ip (x) 

involving G(-) = G (•) + A (x) ■. It holds that 

u{x,t) = -E,(^efo^^-^)^^ij{xt)yo 
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